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Abstract 

I give an introduction to the theory of preasymptotic effects based on the systematic 
OPE/HQET expansion in 1/mn where niQ is the heavy quark mass. The general idea 
is explained in two most instructive examples, with an emphasis on pedagogical aspects. 
Some important results of the last year are reviewed. In discussing the issue of the quark- 
hadron duality, one of the basic ingredients of the theory, I prove that the operator product 
expansion per se is an asymptotic expansion. The behavior of the high order terms in 
this expansion determines the onset of duality and the accuracy of the duality relations. 
The factorial divergence of the high-order terms in OPE implies a sophisticated analytical 
structure in the a s plane, with terms of the type exp[— exp(l/a s )]. 



'Extended version 



1 Introduction 



Although the weak inclusive decays of heavy flavor hadrons are driven by weak interactions 
the theory of pre-asymptotic effects is an important applied branch of QCD. The theory 
has been gradually crystallizing from the eighties, with a decisive breakthrough achieved 
recently, in the last two or three years. I will try to summarize what we know and what 
has to be done in future. 

My talk is preceded by that of Dr. Neubert covering, in part, the same topic. Therefore, 

1 will avoid those aspects that have been already discussed. My task is two-fold. First, I will 
sketch a general picture of the phenomena lying behind the pre-asymptotic corrections. It 
is remarkable that the most essential aspects of the theory that I am going to describe have 
very close parallels with the well-known phenomena from general physics and/or quantum 
mechanics, for instance, the Doppler effect, the Mossbauer effect, etc. It is a pleasure to 
reveal these analogies by translating from the language of the QCD practitioners to that of 
general physics. 

As it often happens, practical needs make us think of deep theoretical questions. One 
of such issues, closely related to the theory of preasymptotic effects, is the onset of duality 
in QCD. When and how it sets in? The answer to these questions is important for under- 
standing the accuracy of numerous predictions of our theory. Some ideas relevant to this 
issue will be discussed in the end of my talk. 

The second task is to present the theory, as it exists now, in a historical perspective. 
The majority of the theorists active in the field are so young that they perceive 10-year old 
results as an ancient history. Therefore it would be in order to remind that both, the ideas 
we rely upon now and the relevant mathematical apparatus, have emerged in the eighties 
being, in turn, a natural continuation of the previous work. 

2 Formulation of the problem 

Heavy flavor hadrons Hq contain a heavy quark Q (practically one can speak of the b 
quark and, to a lesser extent, of the charmed quark) plus a light cloud built from light 
quarks (antiquarks) and gluons. The heavy quark Q experiences a weak transition. The 
nature of this transition is of no concern to us here. It can be a radiative transition, like 
b — ► sj, semileptonic decay like b — ► clv\ or a non-leptonic decay to lighter quarks, e.g. 
b — > cud (Fig. 1). It is assumed that at short distances the amplitude is known from the 
electroweak theory. The task is to calculate the decay rate of the hadron Hq and other 
decay characteristics: the energy spectra, the average invariant mass of the hadronic state 
produced, etc. 

Asymptotically, if the heavy quark mass itiq is infinitely large, the presence of the light 
cloud around Q plays no role. One can treat Q as a free heavy quark decaying into free 
final quarks. This is the famous parton approximation of the pre-historic era, the seventies. 
Having accepted this prescription we, of course, immediately calculate whatever we are 
asked to. Moreover, the result is obviously universal, independent of the structure of the 
light cloud. In particular, all lifetimes are predicted to be the same, 

T(H Q )=T(Q). (1) 

While this simple approach is certainly valid in the asymptotic limit mg = oo, and 
the universality is the feature of this limit, it is totally unsatisfactory for down-to-earth 



1 



physics of the actual heavy quarks. Indeed, a single glance at the lifetime hierarchy in the 
charmed family shows how far we are from the asymptotic regime. The D + meson, the 
most long-living member of the family, has the lifetime exceeding that of by more than 
an order of magnitude, 

t(D+)/t(~° c ) ~ 12. (2) 

Thus, one can hardly speak of pre-asymptotic corrections in this case - the "corrections" 
actually totally determine the decay pattern. In the b quark family they are expected 
to be more modest, but still theoretical control over the inclusive decays is impossible 
without calculating pre-asymptotic effects in a systematic manner. Their dynamical origin 
is quite obvious: interaction of the heavy quark Q and its decay products with the light 
cloud produces a series of terms suppressed by inverse powers of tuq. (Perturbative QCD 
corrections emerging from the hard gluon exchanges and suppressed logarithmically, by 
powers of a s (mg), can be trivially incorporated in the parton model and will not be touched 
upon here.) 

3 Quarks and gluons versus general physics 

Different dynamical mechanisms generate the 1/rriQ corrections. The essence of the phe- 
nomenon can be best explained by a transparent example familiar to everybody from the 
standard courses on quantum mechanics and nuclear physics. Assume we are interested in 
the nucleon (3 decay; the nucleon at hand is not free but is rather bound inside a nucleus. 
The nucleus contains several protons and neutrons glued by the nuclear forces and is sur- 
rounded by an electron shell (Fig. 2). The energy release in the (5 decay, AE, is certainly 
much larger than a typical binding energy of the electrons in the shell. Moreover, we will 
further assume - counterfactually - that AE 3> e nuc i where e nuc i is the nucleon binding 
energy in the nucleus. Such nuclei do not exist in nature, but for our pedagogical purposes 
it is convenient to pretend that they do. 

We start from a free neutron (5 decay, n — > pev e , and then ask what distinguishes it 
from that of a bound nucleon. The answer is well-known: (i) the Pauli interference; (ii) the 
K capture (by proton); and, finally, (iii) the Fermi motion, i.e. the momentum spread of 
the decaying neutron. Let us discuss these phenomena in turn. 

The Pauli exclusion principle does not allow the electron produced to have the same 
momentum as that of the electrons in the shell. Typically the electron momentum is of 
order AE, much larger than that of the bound electrons, and interference does not occur. 
In a "corner" of the phase space, however, the electron produced is soft implying a strong 
interference. This happens rarely; the effect is clearly suppressed by inverse powers of AE. 

The inverse (3 decay due to the K capture is simply impossible for a free nucleon. The 
power suppression here is due to the fact that the amplitude is proportional to Y>(0), the 
bound electron wave function at the origin. 

Finally, the momentum spread of the bound nucleon - a typical bound state effect - 
results in a smearing of the momentum distribution of the produced electron. If the decaying 
neutron has the momentum parallel to that of the electron the latter is accelerated while 
for the antiparallel configuration it is decelerated. As a result, the total decay rate is also 
distorted, although the distortion is parametrically small provided that e nuc i/AE <C 1. 

The description above is qualitative. Quantitative analysis in the case at hand presents 
no difficulties. The standard quantum-mechanical methods and the knowledge of the wave 
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functions (of the electrons in the shell as well as the nucleons in the nucleus) solves the 
problem. For instance, to take into account the Pauli interference we just antisymmetrize 
the total wave function with respect to the coordinates of the electron produced and that 
in the shell. In this way we get the interference term automatically. To include the Fermi 
motion the neutron momentum distribution is convoluted with the amplitude of the free 
decay. 

The pre-asymptotic effects in the inclusive decays of heavy flavors are perfectly similar 
in nature. The decay products of the heavy quark can interfere with the quarks from the 
light cloud. The heavy quark can capture a light one from the cloud. Finally, the "Fermi 
motion" of Q is also there: even if Hq is at rest Q has a non-vanishing spatial momentum 
of order of Aqqd- What is drastically different, however, is the fact that the wave function 
approach is useless in QCD. 

Indeed, the light cloud surrounding Q is a messy dynamical system composed of an 
infinite number of degrees of freedom, quark and gluon, and any attempt to reduce its 
dynamics to quantum mechanics of a finite number of particles is doomed to failure from 
the very beginning. 

Does this mean that no adequate theory can be developed? The question needs no 
answer since the answer has been already spelled out. A systematic approach does exist, 
but one has to pay a price. If in the quantum-mechanical case the effects listed above are 
tractable for arbitrary values of AE - the requirement AE 3> e n ud was actually unnecessary 
- in the heavy flavor decays the theory is based on the expansion in AqcWttiq; practically 
only the first few terms are calculable. Still, even this limited progress in the quantitative 
description of the inclusive decays is of great practical value and leads to a remarkably rich 
phenomenology. 

4 Formalism. A little bit of history 

As usual in practical problems from quantum chromo dynamics the only tool existing today 
for dealing (analytically) with non-perturbative effects is the Operator Product Expansion 
(OPE) [|J. The general idea of OPE was formulated by Wilson even before QCD and is 
very simple. Whatever amplitude is considered all contributions to this amplitude must be 
systematically sorted out into two classes: hard contributions and soft ones. One introduces 
an auxiliary parameter, fj,, (normalization point); if characteristic frequencies in the given 
fluctuation are higher than \x this fluctuation is to be ascribed to the hard part and deter- 
mines the coefficient functions. In QCD the coefficient functions are explicitly calculable 
provided that /i is chosen appropriately, 

a s (n)/ir < 1. 

(Warning: one should not think that the coefficient functions are determined exclusively 
by perturbation theory; in principle, non-perturbative contributions may also be present, 
for instance, those from small-size instantons. What is important is the fact that the 
weak coupling regime is applicable to the coefficient functions. Practically, in the vast 
majority of problems, the coefficient functions are given by perturbation theory to a good 
approximation. This approximate rule is sometimes called the practical version of OPE ||.) 

Soft fluctuations are referred to the matrix elements of local operators appearing in OPE. 
Their contribution is not analytically calculable in the present-day QCD. The corresponding 
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matrix elements can be parametrized, however, or extracted from other sources by using 
symmetry arguments, QCD sum rules or lattice calculations, etc. All dependence on large 
parameters, such as tbq, resides in the coefficient functions. 

The general strategy outlined above has its peculiarities in application to the inclusive 
weak decays. The standard procedure is formulated in the euclidean domain where the 
distinction between the short and large distances (high and low frequencies) causes no 
doubts. Kinematics of the heavy flavor decays is essentially Minkowskian, however; the 
large expansion parameter is the energy release rather than the Euclidean off-shellnes. The 
light quarks, once produced in the Q decay, stay forever - thus, formally, we have to deal with 
the infinite time intervals. Large distance dynamics manifests itself in full in the exclusive 
decays. At the same time, in the inclusive description, after summing over all possible 
channels and certain energy integrations, the results for the coefficient functions depend 
only on the short distance dynamics. This assertion, known as the quark-hadron duality, 
can be justified by means of an analytic continuation in the complex plane in auxiliary 
momenta. The accuracy of the quark-hadron duality is exponential in large momenta 
relevant to the given problem. I will return to this point, the exponential accuracy, later 
on. 

Technically, in order to launch the OPE-based program one has to consider the so-called 
transition operator T related to the modulus squared of the amplitude of interest ||, rather 
than the decay amplitude per se, 



where £w is the short-distance weak lagrangian governing the transition Q — > f under 
consideration. The finite state / can be arbitrary; Fig. 3 displays, as a particular example, 
the graph pertinent to the semileptonic transition Q — ► qlv\ where q stands for the light 
quark. This graph is rather symbolic; one should keep in mind that only the skeleton lines 
are depicted, and all these lines are actually submerged in a soft gluonic medium coming 
from the light cloud in Hq. 

Next, we use the fact that the momentum operator of the heavy quark Q contains a 
large c-number part, 



where is the four- velocity of the heavy hadron Hq, ir^ is the residual momentum operator, 



and G^u is the gluon field strength tensor corresponding to the "background" gluon field 
in the light cloud. The presence of the large mechanical part in the Q quark momentum 
operator ensures that all lines encircled by a dashed line on Fig. 3 can be treated as hard. 
This block, thus, shrinks into a point; it determines the coefficient functions in front of local 
operators built from soft lines. The the set of the soft lines includes the gluon field, light 
quark fields and those describing Q. 

Once the general idea is understood it does not take much effort to calculate T in the 
form of the Wilson expansion, 




(3) 



(4) 





(5) 
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where the operators O n are ordered according to their dimensions (or, under certain circum- 
stances, according to their twists). The coefficients C n are proportional to the corresponding 
powers of 1/mg, modulo calculable logarithms. 

At the next stage we average the transition operator over the hadronic state Hq, i.e. 
find (Hq\T\Hq) . It is at this stage that the internal structure of Hq enters in the analysis. 
The matrix elements {Hq\O u \Hq) parametrize the dynamical properties of the light cloud 
and its interaction with the heavy quark Q. They substitute the wave function of the 
quantum-mechanical treatment, the closest substitute one can possibly have in QCD. 

Finally, at the very last stage, to make contact with measurable quantities we take the 
imaginary part of the matrix element (Hq\T\Hq) . In the example suggested above (Fig. 3) 
this imaginary part is proportional to the total semileptonic decay rate, 

(M HQ r 1 Im{H Q \f\H Q ) = T S L. (6) 

If the lepton energy is not integrated over but is held fixed the imaginary part of the 
transition operator is proportional to the energy spectrum, etc. 

Analysis of the weak inclusive decays along these lines has been suggested in the mid- 
eighties when some operators in the expansion (5) relevant to the total lifetimes in the 
families of charm and beauty have been considered ||. Many elements of what is now 
called the Heavy Quark Effective Theory (HQET) £§ § can be found in these works. 
Further progress has been deferred for a few years by an observation of certain difficulties 
in constructing OPE in the Minkowski domain || . The difficulties turned out to be illusory 
and are explained by the fact that contributions discussed in Ref. Q were not fully inclusive. 
Once this has been realized Q the last obstacle disappeared and a breakthrough in the 
theory of pre-asymptotic effects became inevitable. First systematic results - the leading 
1 / ttiq corrections to the parton- model predictions - have been obtained shortly after |8], ||] . 
While the works of the eighties passed essentially unnoticed fusion of similar ideas, which 
took place, independently, in the very beginning of the nineties, resulted in HQET in its 
present formulation j|, || . HQET provided a simple and convenient language for discussing 
the l/m,Q expansions, in particular, in connection with the semileptonic decays. General 
aspects of the HQET/OPE approach to the semileptonic inclusive decays of heavy flavors 
have been discussed in Ref. |l(J . The authors of this important work were not familiar with 
the previous development. 

The current state of the theory of pre-asymptotic effects is characterized by a rapid 
expansion in the practical spheres. The 1/rriQ, 1/rriQ and even 1/r/ig corrections are found 
for probably all conceivable quantities measurable in the inclusive decays of and H c . 
An incomplete list of works published in the last year and devoted to this subject includes 
a dozen of papers [11-22]. What remains to be done to complete the picture? The initial 
excitement, quite natural in time of the explosion of the theoretical predictions based on 
first principles, should give place to an every-day routine of extracting reliable and model- 
independent numbers from experimental data. This aspect - relation to experiment - will 
not be touched upon in my talk. Instead, I will demonstrate in more detail how the theory 
works in two problems which seem to be most instructive. 

5 Total semileptonic widths. CGG/BUV theorem 

The inclusive semileptonic decay rate is probably the simplest and the cleanest example. 
Although the parton-model analysis (including the first radiative QCD correction) is known 
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for more than a decade |23| the leading non-perturbative correction was calculated only in 
1992, a pilot project ]|] which opened the modern stage and was supplemented later by 
calculations of the energy spectrum and other distributions. 

If one neglects the hard gluon exchanges, the transition operator is given by the diagram 
depicted on Fig. 4. The weak Lagrangian responsible for the semileptonic decay has the 
form 

C = ^V Qq {qT li Q)(ir il u), (7) 

where I is a charged lepton, electron for definiteness, and we assume for simplicity that the 
final quark q is massless (T^ = 7^(1 + 75))- It is understood that the quark lines on Fig. 4 
describe propagation in a soft gluon medium. The leptons are, of course, decoupled from 
gluons, and their lines are always treated as free. 

The operator product expansion for the transition operator T starts from the lowest- 
dimension operator, QQ in the case at hand. In calculating the coefficient of this operator 
one disregards any couplings of the quark lines to the soft gluons from the cloud. The full 
perturbation theory resides in this coefficient, but - as I will explain shortly - it is wrong to 
say that the QQ term in OPE just reproduces the parton-model prediction for Fsl] actually 
it gives more. 

There are no gauge and Lorentz invariant operators of dimension four (and only 
Lorentz and gauge invariant operators can appear in the expansion for Tsl), an d the first 
subleading term in OPE is due to the dimension 5 operator 

O g = Q'-o^G^Q. (8) 

This is the only Lorentz scalar dimension-5 operator bilinear in Q. 

The corresponding coefficient Cg can be calculated in a variety of ways, of course; the 
easiest one is based on the background field technique and the Fock-Schwinger gauge [Q. 
In this gauge the line corresponding to the final quark q remains free (provided that we limit 
ourselves to the operator Oq), and the only reason why the dimension-5 operator appears 
at all is the fact that we use the exact Dirac equations of motion with respect to the heavy 
quark lines Q and Q. 

For those who are not familiar with the Fock-Schwinger gauge technique let me add 
a remark elucidating the above statement. Explicit calculation of the diagram of Fig. 4, 
where all quark lines are considered in the background gluon field (with the Fock-Schwinger 
gauge condition), yields the operator Q(0) p / (p 2 ) 2 Q{0) where p = id. Now, we would like to 
rewrite p in terms of the covariant derivative V = p + A and the gluon field strength tensor 
- which is perfectly trivial in the Fock-Schwinger gauge - and then take advantage of the 
fact that IjPQ = itiqQ. It is not difficult to see that in the linear in G^ u approximation one 
can substitute 

Q(0) y(p 2 ) 2 Q(o) -» Q(0) fCP 2 ) 2 Q(o). 

Supplementing this expression by the equalities 

V 2 = f 2 - % -oG and fQ = m Q Q 

we conclude that in the Fock-Schwinger gauge the operator Q fl{p 2 ) Q reduces to itiqQQ — 
2tjiqOg- In other words, calculating the diagram of Fig. 4 with free fermions we get, for 
free, both coefficients - the coefficient in front of QQ and that in front of Og^ 
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I bring my apologies for this brief technical digression intended to demonstrate the 
beauty of the Fock-Schwinger technique which, unfortunately, is still unfamiliar to many 
heavy quark practitioners. Now it is time to present the final prediction || for the inclusive 
semileptonic decay rate, 

Glm 5 n \Vna\ 2 

1 2M ffQ m\ 2M Hq 1 « 'J 

Perturbative hard gluon corrections are omitted here. 

This is almost the end of the story. To complete the analysis one needs to know the 
matrix elements of the operators QQ and Oq over the hadronic state Hq of interest. These 
matrix elements, clearly, depend on the structure of the light cloud which is beyond our 
computational abilities at present. Fortunately, they are known independently! 

Indeed, Og is the lowest-dimension operator responsible for the spin splittings - the 
mass difference between the vector and pseudoscalar mesons (e.g. B* and B). Thus, 

^^ = »G = \{Mt-M*) (10) 

where M* and M are the B* and B masses, respectively. The right-hand side of eq. (10) 
is well measured experimentally. As for (Hq\QQ\Hq), let us compare the operator QQ 
with the current Q^oQ whose matrix element over any Q flavored hadron is unity, of course 
(in the Hq rest frame). We are lucky again. The difference QQ — Q^/oQ is due to the 
lower components of the bispinors, and for heavy quarks these are reducible to the upper 
components by means of the equations of motion expanded in 1/rriQ. In this way we arrive 
at§ 

QQ = Q l0 Q + J- Y O G --^ Y O 7T + O(m Q 4 ), (11) 

ZITIq ZUIq 

where 

= Q^Q 

has the meaning of the average spatial momentum squared of the heavy quark Q inside Hq. 

As you probably remember, I told you that all perturbation theory is hidden in the 
coefficient of the operator QQ. The parton- model result is recovered if one approximates 
the matrix element (2M^) _1 (Hq\QQ\Hq) by unity. Now we see that this is actually the 
zeroth order approximation; the very same operator QQ gives rise to \jm 2 Q and higher-order 
terms, for instance, 

— !— (B\bb\B) = l + -^--^E_ + ... (12) 
2Mb 2ml 2m b 

where ^ parametrizes the matrix element of O^, 



^ = (2M HQ r 1 (H Q \Q7f 2 Q\H Q ). 
Assembling all these results together we finally find [||] 

4m I 



T { B^X M - i927r3 
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(13) 



where the expression in the first square brackets is due to the 1 /m% correction in the matrix 
element of bb, the term in the second square bracket comes from the operator Oq in OPE 
for Vshi and the dots denote higher-order terms in l/m^. 

The //^ part in the first square bracket has a very transparent physical interpretation. 
This is nothing else than the quadratic Doppler effect reflecting the time dilation for the 
heavy quark decaying in flight. Indeed, in the rest frame of the decaying hadron the quark 
b has a spatial momentum prolonging its lifetime (diminishing T). The effect linear in this 
momentum drops out because of averaging over angles, while the quadratic effect survives. 
It is quite obvious, from the first glance, that the coefficient of ^ in eq. (13) is exactly as 
it is expected to be in this picture. 

A remarkable fact which I have especially emphasized in eq. (13) is the absence of the 
correction linear in 1/rriQ in the total semileptonic width. This practically important theo- 
rem (CGG/BUV theorem) is derived from two components - the absence of the dimension-4 
operators in OPE and the absence of 1/rriQ corrections in the matrix elements of the op- 
erator QQ. The first circumstance has been noted in Ref. 10. The authors of this work 
believed, however, that linear in 1/mg terms may appear in (2M) _1 (Hq\QQ\Hq). The 
non-renormalization theorem for (2M) _1 (Hq\QQ\Hq) at the 1/rriQ level has been estab- 
lished in Ref. 8 completing the proof of the theorem. 



6 The impact of the heavy quark motion in Hq 

Even if the Q flavored hadron is nailed at the origin so that its velocity vanishes the heavy 
quark Q moves inside the light cloud, its momentum being of order Aqcd- This is the QCD 
analog of the Fermi motion of the nucleons in the nuclei. It is quite clear that this motion 
affects the decay spectra. Say, if the "primordial" heavy quark momentum is parallel to that 
of the lepton pair, the pair gets boosted. A distortion of spectra may lead to corrections in 
the decay rate. It is equally clear that this effect is pre-asymptotic (suppressed by inverse 
powers of toq): while typical energies of the decay products are of order mg a shift due to 
the heavy quark motion is of order Aqcd- 

A few models engineered to describe this effect at a purely phenomenological level exist 
in the literature for years (e.g. Ref. 25). Now one can finally approach the issue scientifically, 
backed up by a solid QCD-based analysis. The theory of this phenomenon is very elegant. I 
would like to demonstrate basic points in a simplest example, Q —* qj transition. Moreover, 
to avoid inessential technicalities (which can be readily worked out, though) I will neglect 
the quark and photon spins, assuming that both fields, Q and q are spinless. As for the 
mass of the final quark m q it will be treated as a free parameter which can vary from zero 
almost up to uiq. (For our approach to be valid we still need that Am = m,Q — m q 3> Aqcd 
although the mass difference may be small compared to the quark masses. As a matter of 
fact, this small velocity (or SV) limit [26] is very interesting and instructive, and I will 



return to its discussion later on.) 

Thus, for pedagogical purposes we will consider, following Ref. 17, a toy model, 

£ = hQ(f>q + h.c. (14) 
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where (ft is a toy "photon" field and h is a coupling constant. To begin with, it is convenient 
to put the final quark mass to zero. At the level of the free quark decay the photon energy 
is fixed by the two-body kinematics, = mq/2. In other words, the photon energy 
spectrum is a monochromatic line at E$ = hiq/2 (Fig. 5). On the other hand, in the actual 
hadronic decays the kinematical boundary of the spectrum lies at Mh q /2. Moreover, due to 
multiparticle final states (which are, of course, present at the level of the hadronic decays) 
the line will be smeared. In particular, the window - a gap between mq/2 and Mh q /2 - will 
be closed (Fig. 6). There are two mechanisms smearing the monochromatic line of Fig. 5. 
The first is purely perturbative: the final quark q can shake off a hard gluon, thus leading to 
the three-body kinematics. The second mechanism is due to the "primordial" motion of the 
heavy quark Q inside Hq and is non-perturbative. Only the second mechanism will be of 
interest for us here. The theory of the line shape in QCD resembles that of the Mossbauer 
effect. 

Thus, the question to be answered is: "how can one translate an intuitive picture of the Q 
quark primordial motion inside Hq in QCD-based theoretical predictions, not just plausible 
models?" The OPE /HQET approach is applicable here in full, with a few additional - rather 
transparent - technical elements. 

The basic object of consideration is again the transition operator (3) (with Cw substi- 
tuted by Since we are interested in the energy spectrum the "photon" momentum q 
must be fixed. 

Then in the approximation of Fig. 7 (no hard gluon exchanges) the transition operator 
takes the form 

f = - f d 4 x(x\Q 1 Q|0), (15) 

where the background field technique is implied and (Po)n = mQV^. I remind the toy model 
with no spins for the quark and photon fields is taken for illustration; bellow a convenient 
shorthand notation 

k = P - q = niQV - q 

will be used. 

If T is calculated in the operator product expansion the problem is solved since 

rfV h 2 E 

liE = ^M^^ H ^ < 16 > 

To construct OPE we observe that, as in the previous problem of the total semilep- 
tonic widths, the momentum operator tt corresponding to the residual motion of the heavy 
quark is ~ h-QCD and the expansion in ir/k is possible. Unlike the problem of the total 
widths, however, now in the interesting energy range k 2 is anomalously small, the expansion 
parameter is of order unity, and an infinite series of terms has to be summed up. 

To elucidate this statement let us examine different terms in the denominator of the 
propagator, 

k 2 + 2-nk + -K 2 . 

The effect of the primordial motion plays a key role in the window and in the adjacent 
region below it, 

\E — (l/2)rriQ\ ~ A (17) 
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where A = Mh q — tuq. It is quite trivial to find that in this domain 

ko ~ \k\ ~ mg/2, k 2 ~ tuqA; 

in particular, at the kinematical boundary (for the maximal value of the "photon" energy) 
fco = Mq/2 and k 2 = —tuqA. Hence, inside the line 

k 2 ~ kir ^> 7r 2 . 

In other words, when one expands the propagator of the final quark in eq. (15) in tt, in 
the leading approximation all terms (2kn/k 2 ) n must be taken into account while terms 
containing tt 2 can be omitted. The first subleading correction would contain one 7r 2 and 
arbitrary number of 2/c7r's, etc. For those who still remember the theory of deep inelastic 
scattering (DIS) of the late sixties and early seventies this may sound suspiciously familiar. 
Yes, you are right; the theory of the line shape in the radiative transition of the heavy quark 
into a massless one can be viewed as a version of the standard theory of DIS. 

Thus, in this problem it is twist of the operators in OPE, not their dimension, that 
counts. Keeping only those terms in the expansion that do not vanish in the limit niQ — * oo 
(analogs of the twist- 2 operators in DIS) we get the following series for the transition 
operator JP7| ] 

1 00 f 2 \ n 

r = fc w ...fc^(Q7r M1 ...7T Aln Q- traces). (18) 

Traces are subtracted by hand since they are irrelevant anyway; their contribution is sup- 
pressed as k 2 ir 2 /{kir) 2 ~ A/rriQ to a positive power. Another way to make the same 
statement is to say that in eq. (18) the four-vector k can be considered as light-like, k 2 = 0. 

After the transition operator is built the next step is averaging of T over the hadronic 
state Hq . Using only the general arguments of the Lorentz covariance one can write 

{H Q \QTT fll ...TT fln Q - traces|iiQ) = a n A (v^.-.v^ - traces) (19) 

where a n are constants parametrizing the matrix elements. Their physical meaning will 
become clear momentarily. Right now it is worth noting that the term with n = 1 drops 
out (ai = 0). Indeed, {Hq\QttQ\Hq) is obviously zero for spinless Hq while ttq through the 
equation of motion reduces to -k 2 j(2mQ) and is of the next order in 1/mg. Disappearance 
of Qtt^Q means that there is a gap in dimensions of the relevant operators. 
Let us write a n 's as moments of some function F(x), 

a n = J dxx n F{x). (20) 

Then, F(x) is nothing else than the primordial line-shape function! (That is to say, F{x) 
determines the shape of the line before it is deformed by hard gluon radiation; this latter 
deformation is controllable in perturbative QCD). The variable x is related to the photon 
energy, 

x = (2/A)[E - (l/2)m Q ]. 

If this interpretation is accepted - and I will prove that it is correct - it immediately implies 
that (i) F(x) > 0, (ii) the upper limit of integration in eq. (20) is 1, (iii) F{x) exponentially 
falls off at negative values of x so that practically the integration domain in eq. (20) is 
limited from below at — xq where xq is a positive number of order unity. 
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To see that the above statement is indeed valid we substitute eqs. (19), (20) in T, 

mfm - -i e/ *fm »" • < 21 > 

and sum up the series. The iO regularization will prompt us how to take the imaginary part 
at the very end. In this way we arrive at |l7j 

where the variable x is defined above and To is the total decay width in the parton approx- 
imation. Corrections to eq. (22) are of order A/mg. 

Thus, we succeeded in getting the desired smearing: the monochromatic line of the 
parton approximation is replaced by a finite size line whose width is of order A. The 
pre-asymptotic effect we deal with is linear in Agco/mg- 

At this point you might ask me how this could possibly happen. I told you a few minutes 
ago that there is a gap in dimensions of operators in OPE - no operators of dimension 4 
exist - and the correction to QQ is also quadratic in 1/mn (the CGG/BUV theorem). 
No miracles - the occurrence of the effect linear in Aqcd/it^q became possible due to the 
summation of the infinite series in eq. (18); no individual term in this series gives rise to 
^QCD/mQ. 

To avoid misunderstanding it is worth explicitly stating that the primordial line function 
is not calculated; rather F(x) is related to the light-cone distribution function of the heavy 
quark inside Hq, vis. (Q(mr) n Q), n 2 = 0, or more explicitely 

F(x) oc J dte lxiK {H Q \Q{x = 0)e~ i So nA ^ dr Q(x^ = n^)\H Q ), (23) 

where n is a light-like vector 

n (t = (1,0,0,1). 

Moreover, this primordial function is not the one that will be eventually measured from 
dT/dE; the actual measured line shape will be essentially deformed by radiation of the 
hard gluons. I will say a few words about this effect later on. 

Although we can not calculate F(x) still our analysis is useful in many respects. First of 
all, a few first moments of F{x) are known, as well is its qualitative behavior. Second, the 
QCD-based approach tells us that the "scientific" distribution function is one-dimensional, 
while in a naive model one would rather introduce a three-dimensional distribution in the Q 
spatial momentum. The reason is obvious: the momentum operator 7r M does not commute 
with itself, [vr^Tr^] = iG^ u , while nir commutes. Finally, the logarithmic evolution of F(x) 
is calculable in perturbative QCD. 

For real quarks and photons, with spins switched on, calculation of the line shape is 
carried out in Refs. 17,18. 

7 Varying the mass of the final quark 

So far I was discussing the transition into a massless final quark. If we looked at the line 
shape presented on Fig. 5 more attentively, through a microscope, we would notice that a 
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smooth curve is obtained as a result of adding up many channels, specific decay modes. A 
typical interval in E that contains already enough channels to yield a smooth curve after 

2 

summation is ~ A /mq. In other words, in terms of the photon energy the duality interval 
is ~ A 2 /rriQ. Roughly one can say that the spectrum of Fig. 5 covers altogether mq/A 
resonance states poduced in the Hq decays and composed of q plus the spectator (I keep 
in mind here that the final hadronic state is produced through decays of highly excited 
resonances, as in the multicolor QCD). These states span the window between mq/2 and 
Mhq and the adjacent domain to the left of the maximum at E = mq/2. 

It is very interesting to trace what happens with the line shape, the duality interval and 
the primordial distribution function as the final quark mass m q increases | jT^| . 

The most obvious kinematical effect is the fact that the window shrinks. When we 



eventually come to the SV limit [26| 



Aqod < Am = mq - m q < mq A 

it shrinks to zero. Indeed, in this limit the photon energy in the two-body quark de- 
cay, Am(l + Am(2mg) _1 ) differs from the maximal photon energy in the hadronic decay, 
AM(1 + AM(2Mq) _1 ), only by a tiny amount inversely proportional to mq (Am and AM 
stand for the quark and meson mass differences, respectively). Simultaneously with the 
shrinkage of the window the peak becomes more asymmetric and develops a two-component 
structure (Fig. 8). The dominant component of the peak, on its right-hand side, becomes 
narrower and eventually collapses into a delta function in the SV limit. A shoulder de- 
velops on the left-hand side; the duality interval becomes larger and the number of the 
hadronic states populating the line becomes smaller. When we approach the SV limit the 
delta function component will consist of exactly one sate - the lowest-lying q containing 
meson (the elastic component) while the shoulder (the inelastic component) will include 
several excitations and will stretch down to E = E max — several units xA. The light-cone 
distribution function will evolve and will continuously pass into the temporal distribution 
function determining the form of the spectral line in the shoulder (Fig. 9). 

This rather sophisticated picture, hardly reproducible in the naive models, stems from 
the same analysis [|l7). The transition operator T now has the form 



_y 

ml — k 2 ^ n 



I ^ ^ / 2mqi r + 7r 




k- - 1 * <24> 

where q is the momentum. Notice that 2mqiTQ + tt 2 acting on Q yields zero (the equation 
of motion) and in the SV limit q must be treated as a small parameter, 

Qo/mq = E/mq = v < 1; 

v is the spatial velocity of the heavy quark produced. Although v is small the inclusive 
description we develop is still valid provided that Am 3> A-qcD- 

In the zeroth order in q the only term surviving in the sum (24) is that with n = 0, 
and we are left with the single pole, the elastic contribution depicted on Fig. 9. This is 
the extreme realization of the quark-hadron duality. The inclusive width is fully saturated 
by a single elastic peak. What might seem to be a miracle at first sight has a symmetry 
explanation - the phenomenon is explained by the heavy quark symmetry which has been 
noted in this particular context in Refs. 26, 27. (The heavy quark symmetry is also called 
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the Isgur-Wise symmetry [£8|]. These authors generalized the observation of Refs. 26, 27 to 
arbitrary values of v.) On the other hand, the fact that the parton-model monochromatic 
line is a survivor of hadronization is akin to the Mossbauer effect. 

If terms 0(v 2 ) are switched on the transition operator acquires an additional part, 

From this expression it is obvious that the shape of the v 2 shoulder is given [17] by the 
temporal distribution function G(x) whose moments are introduced through the matrix 
elements 

(Hq I QiTiiTQiriQ | Hq) = A n+2 J dxx n G(x). (26) 
Alternatively, G(x) can be written as a Wilson line along the time direction, 

G(x) oc J dte ix{K (H Q \Q(t = 0, x = O^e"* Ao(T)< %Q(i, x = 0)\H Q ). (27) 

Intuitively it is quite clear why the light-cone distribution function gives place to the 
temporal one in the SV limit. Indeed, if the massless final quark propagates along the 
light-cone, for Am <C m the quark q is at rest in the rest frame of Q, i.e. propagates only 
in time. 

In terms of G{x) our prediction for the line shape following from eq. (25) takes the form 



dT 

dE* 



S(x)+ 



3 \x z x I 

where x = (E — E max )/A. The v 2 corrections affect both, the elastic peak (they reduce 
the height of the peak) and the shoulder (they create the shoulder). The total decay rate 
stays intact, however: the suppression of the elastic peak is compensated by the integral 
over the inelastic contributions in the shoulder. This is the Bjorken sum rule p9| , pCp . It is 
important that we do not have to guess or make ad hoc assumptions - a situation typical 
for model-building - the theory itself tells us what distribution function enters in this or 
that case and in what particular way. 



8 Inclusive semileptonic decayst 161719 ' 21 ] 

The very same distribution functions that determine the line shape in the radiative tran- 
sitions appear in the problem of the spectra in the semileptonic decays. In particular, in 
b — ► ulv we deal with F(x). Their analysis is absolutely crucial if one addresses the behavior 
in the so-called end-point region (i.e. the region inaccessible in the parton model and close 
to it, an analog of the window in the radiative transitions). 

Apart from trivial kinematical modifications - for instance, the occurrence of several 
structure functions - the only change is the expression for the variable x. More exactly, as 
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was shown above, for massless final quark in both cases, radiative transitions and semilep- 
tonic decays, 

x = -k 2 /(2Akv), 

but in the latter case q 2 > while in the radiative transitions q 2 = 0, and this seemingly 
insignificant difference is practically important. Indeed, if q 2 = the notion of scaling 
makes no sense - kinematically there is only one variable, the photon energy E, and x oc 
(E — E max)- On the other hand, in the inclusive decays the structure functions depend, 
generally speaking, on two independent variables, k 2 and kv (or q 2 and qo). The fact that 
the structure functions in B — > X u lv are actually functions of a single combination x is a 
very strong statement fully equivalent to the Bjorken scaling in deep inelastic scattering. 
Guesses about a scaling behavior in the inclusive semileptonic decays are known in the 



literature [31]. This is the first time ever we are able to say for sure what sort of scaling 
takes place, where it is expected to hold and where and how it will be violated. 

In the most naive parton model, when gluons are not considered at all, k 2 = 0. The 
primordial heavy quark motion smears k 2 so that actually 

k 2 ~Ak . (29) 

I will sometimes refer to this domain as to "primordial". In this domain the structure 
functions depend on 

x = -jTV/2/co = -A~V/(fc + |jfe|) = -A _1 (/c - (30) 

where in the denominator the difference between ko and \k\ is neglected which is perfectly 
legitimate in the primordial domain (29). In the rest of the phase space this substitution is 
wrong, of course, but there the above scaling is not going to be valid anyway, so there is no 
need to bother. From eq. (30) we see that in the primordial domain the structure functions 
depend on a single light-cone combination 

Qo + \q\, 

rather than on qo and q 2 separately. 

The primordial structure functions fall off - presumably exponentially - outside the 
primordial domain. The hard gluon emissions will populate the phase space outside this 
domain creating long logarithmic tails. The primordial part is buried under these tails. 
Therefore, outside the primordial domain one can not expect that the structure functions 
depend on the single combination qo + \q\. 

Perturbative corrections as well as those due to higher twist operators violate scaling in 
the primordial domain as well. If the latter are generically small the perturbative corrections 
in the end-point region, due to the Sudakov double log suppression, are expected to be large, 
so that practically in the measured structure functions the scaling feature will be implicit. 
Needless to say that perturbative violations of scaling are calculable. 

I will not go into further details which are certainly important if one addresses the 
problem of extraction of V u b from experimental data. Some of them are discussed in the 



current literature [17, |19[, others still have to be worked out. Applications of the theory 
to data analysis is a separate topic going beyond the scope of the present talk. The only 
obvious remark on the visibility of the primordial motion in experimental data seems to 
be in order: the electron energy spectrum in B —* X u ev is probably the worst place for 
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studying F(x). Indeed, here the primordial distribution function enters in an integrated 
form even in the absence of the hard gluon corrections. A smearing of F(x) is unavoidable 
if we integrate over the neutrino variables, and then the effect of F{x) is diluted, and is 
seen to a much lesser extent than in the radiative transitions. In the inclusive semileptonic 
decays a much better place to search is the double differential distribution in q$ and q 2 . 



9 civ 

The formalism described above is fully applicable [17, El], in principle, to the semileptonic 



inclusive transitions b — ► civ where the final quark c can be also treated as heavy, although at 
the same time , m 2 <C m\. The ratio m 2 /m 2 « 0.07 is a small parameter while m 2 /{Km b ) ~ 
1. Technically, the situation becomes much more sophisticated, however, since the type of 
the distribution function reflecting the primordial motion of the b quark inside the B meson 
will depend now on the value of q 2 where q is the total momentum of the lepton pair. 

If q 2 is small, q 2 < m 2 , one recovers [21] the same light-cone function F(x) as in the 
transition b — ► ulv or b — > 57. Modifications are marginal. First, some extra terms explicitly 
proportional to m c /m.b are generated in the structure functions due to the fact that y + m c 
replaces f in the numerator of the quark Green function. Moreover, if in the b —* ulv 
transitions the scaling variable in the end-point domain is 

x = A _1 (g + \q\ ~ m b ), (31) 

in the b — > civ transition it is shifted by a constant term of order 1, 

2 

1 Tfl 

x = A ( qo + \q\- mb ) + =^. (32) 

Finally, subleading (higher-twist) terms in b — > ulv generically give rise to corrections 
of order k 2 /feg ~ A/mj,. In the beauty-to-charm decays these corrections are 0({hmi + 
m 2 c )/m 2 b ). 



As q increases the corrections to the description [21| based on the light-cone function 



grow and eventually blow up when q 2 approaches its maximal value, 

(Lax = (Mb - M D ) 2 , 
since here ko ~ m c , not m&. Practically this happens at 

fq 2 >M B -2M D . (33) 

The domain (33) corresponds to the SV regime which I have already discussed, with fasci- 
nation, in the toy example above. In the SV limit the velocity of H c produced is small, 

2 ^ (M B -M D ) 2 -q 2 
]V] ~ M B M D 

and the light-cone distribution function becomes irrelevant. The primordial motion of the 



b quark is described [17] in this regime by the temporal distribution function G(x) where 



t— 1 



A (qo — Am) 
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Am is the quark mass difference coinciding, to the leading order with Mb — Mb- The 
"corner" of the phase space where we find ourselves close to the SV limit is responsible for 
roughly 1/2 of all events in B — > X c lu. 

Thus, changing q 2 from zero to q^ax results in an evolution of the distribution function 
appearing in theoretical formulae for dT(B —* X c li>), from light-cone to temporal, through 
a series of intermediate distributions. The physical reason for this evolution is quite clear 
- what distribution function is actually measured depends on the parton-model velocity of 
the quark produced in the b decay. 

In calculating the lepton energy spectrum, dT/dEi, one needs to integrate over a range 
of q 2 . In doing so, we inevitably smear all these evolutionary stages. How exactly the effect 
of the primordial motion in dT(B — > X c li>) / 'dE\ looks like - nobody knows at present. An 
attempt to express this effect in terms of the light-cone function F{x) has been recently 
reported One should expect that, given the actual values of the quark masses, the 

expression for dT(B — > X c lv) / dEi obtained in this work will have corrections 0(1). 

10 The Roman distribution function 

The moments of the light-cone distribution function F{x) are related to the matrix ele- 
ments of the operators (19). This infinite set of matrix elements codes information about 
the bound-state structure of Hq. Needless to say that we are not able to calculate this 
infinite set in the present-day QCD; our possibilities are limited to a few first moments. 
Since the desire to get an idea of the primordial distribution function, to visualize F(x), is 
naturally strong it seems reasonable to adopt a pragmatic attitude. Let us try to conjec- 
ture a qualitatively acceptable function, whose properties are compatible with the expected 
features of the primordial distribution function. This conjectured function will depend on 
several free parameters to be fixed from the moments we know. 

No extensive guesswork is needed. Surprisingly, a phenomenological model existing for 
more than a decade will provide us with a quite adequate input for F(x). I mean the 
so-called AC 2 M 2 model |25j for the "Fermi" motion of the heavy quark inside Hq. 

The physical essence of the model is very simple. It assumes that, say, the B meson 
is a compound system consisting of the heavy quark b and a spectator with a fixed mass 
m sp (one of two free parameters of the model). The spectator component has a momentum 
distribution 3>(|p|) where pis a three-dimensional momentum of the spectator. The energy 
of the b quark is Ef, = Mb — (\p\ 2 + m 2 ^) 1 / 2 . The mass of the b quark then has to be equal 
to Eh through order p. Thus, we arrive at the notion of the floating b quark mass 

m f b ^M B -(\p\ 2 + m 2 sp )^ 2 , (34) 

technically a most important ingredient of the AC 2 M 2 model. The b quark with momentum 
—p decays into uev in flight; after averaging over 




we obtain all relevant decay spectra (here pp, a Fermi momentum, is a second free parameter 
of the model). 

This model is routinely used for years in analyzing the experimental data on the inclusive 
semileptonic decays of B 7 s. Given all its naivete it is quite remarkable that it still reproduces 



16 



essential features of the QCD-based theory of pre-asymptotic corrections. First of all, it 



has been demonstrated [32, 33] that the corrections of the first order in 1/mb to the total 
decay rates are absent in this model, the CGG/BUV theorem. Second, in the primordial 
domain the correct scaling is recovered: the structure functions do indeed depend only on 
the combination (30) or (31). 



It is not difficult to extract [33] the light-cone distribution function stemming from the 
AC 2 M 2 model, 



, . 1 A 1 
FRomix) = -= — exp < -- 
y 7T pp 4 



PF p _ A. 
A 1 — x pf 



(35) 



where Rom stands for Roman, and 

P = (m sp /p F ) 2 . 
The vanishing of the first moment implies 

- = ^ 2 K l( f), (36) 

PF V 71 " 2 

where K\ is the McDonald function. The Roman function satisfies the requirements one 
would expect from the primordial distribution function in the inclusive decays into massless 
quarks. First, it vanishes when one approaches x = 1 from below, so that no spectrum is 
generated beyond the physical kinematical boundary. Moreover, at the negative values of 
x it falls off exponentially once |x| 3> 1. 

The Roman function depends on two parameters to be fitted, m sp and pf- Usually they 
are fitted from experimental data. Previously, when the theory of pre-asymptotic effects 
did not exist, there was no clear idea how to fit the data properly. New understanding calls 
for a new analysis, currently under way. Meanwhile it is instructive to see what can be said 
theoretically. 

It is not difficult to show that the AC 2 M 2 model constrains the second moment from 
above, (x 2 ) < (it — 2)/2. The largest possible value, (x 2 ) ~ 0.6 is achieved if p — ► 0. On 
the other hand, from the QCD sum rules one expects |34| (x 2 ) ~ 0.5 to 1. Therefore, if 
the AC 2 M 2 ansatz has chances to survive the value of p must be small. In other words, 
m 2 p <C Pp, the spectator must be relativistic. In this case Fn om is rather broad - its width 
is of order A - and very asymmetric, see Fig. 10. If p is indeed small, pp ~ (^/^^/2)A. 

It may well happen that the Roman ansatz, being qualitatively reasonable in the b — > u 
transition, will fail to quantitatively reproduce fine structure. With two free parameters it 
may turn out to be too restrictive. Then we will have to invent another clever ansatz, with 
three fit parameters. 

I hasten to add that the AC 2 M 2 model does not reproduce correctly the \jm 2 Q and 
higher effects, and it should not - it was not engineered for that purpose. What may 
be even more important, it gives no hints on the evolution of the light-cone distribution 
function towards the temporal one with increasing m q . Therefore, the description based on 
this model is generally speaking inapplicable in the b —* c transition; especially detrimental 
it becomes in the B — > X c ev decays at large q 2 , when we approach the SV limit. As I have 
mentioned a few minutes ago a large part of the phase space belongs to the SV regime. 
The strongest indication is the fact that two elastic modes, B — > Deu and B — > D*ev, yield 
more than 60 % of the total probability. In the SV regime the dominance of the elastic 
modes is absolute. 
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11 A few words about hard gluons 



So far non-perturbative effects were my prime concern. This is the heart of the theory, its 
non-trivial part. To make contact with experiment, however, it is necessary to include a 
dynamically rather trivial effect due to the hard gluons, in the bremsstrahlung and in the 
loops. The bremsstrahlung corrections are most important in the case of the light final 
quarks and are moderate if the final quark is relatively heavy. Therefore, to make my 
point, I will consider, as a most instructive example [17], the radiative transition b — » S7 



(B — ► Xs^y at the hadron level). As you surely remember, in the parton approximation 
the photon spectrum is a monochromatic line at E = mi/2. The primordial motion of the 
heavy quark, a non-perturbative effect, smears this line and generates a spectrum in the 
window and in the adjacent region below the window. The width of the line becomes of 
order A. Outside the primordial domain whose size is ~ A the primordial motion has a 
negligible impact on the photon spectrum. 

Even if the primordial motion is switched off, the line will be smeared anyway due 
the hard gluon emission. A bremsstrahlung gluon can carry away a finite fraction of the 
accessible momentum, ~ This contribution produces events with the invariant mass 

of the hadronic state of order nib/2; the corresponding photon energy spectrum is 0(a s ). 
The effect is calculable as a one-loop correction in the transition operator. 

Apart from this trivial radiative tail the hard gluon emission modifies the peak itself, 
much in the same way as the electromagnetic radiative corrections modify the J /if) peak 
in e + e~ annihilation. There the natural width of the J/ip meson is negligibly small, and 
the observed shape of the peak is determined by the radiative smearing due to the photon 
bremsstrahlung. 

The gluon bremsstrahlung is characterized by large logarithms in the end-point domain; 
summation of these logarithms results in the Sudakov exponent. The Sudakov effect destroys 
the monochromatic line of the parton model. (It does not populate the window, though; it 
clearly remains empty in perturbation theory). 

The peak at E = m^/2 becomes less singular but still a singularity persists. The total 
probability remains the same as in the parton model but, if in the parton model it is all 
saturated in the peak, now a part of it is pumped into the adjacent tail. The area under 
the distorted peak turns out to be ~ (A/mb) e ° where eo is a positive but relatively small 
number, eo ~ 0.3 |L7|], 

For some reasons which I do not understand the problem of the Sudakov suppression 
in the end-point domain is perceived as something intractable: some authors say that the 
running coupling constant controlling these corrections blows up and no reliable calculations 
are possible. In principle, these questions have been studied in similar problems in the late 
seventies. To feel the running nature of a s it is necessary to go beyond the double-log level. 
The answer is well-known at least at the first subleading level (with one logarithm lost): 
the coupling constant entering the Sudakov exponent is the running coupling normalized 
at k 2 where k 2 is the invariant mass squared of the hadronic state in the corresponding jet. 
In our case of the b — ► 57 transition k 2 = (mju — q) 2 . 

If so, the transformation of the parton-model delta function is easy to find. The Sudakov 
logarithms generate the following spectrum at the end-point: 

dT/dE cx -. — -. - T - 0(— - E) (37) 

' [1 - {2E/m b )} 1 -^ y 2 ' K 1 
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where eo = 8/(36) ~ 0.3 (b is the first coefficient in the Gell-Mann-Low function). This 
expression neglects subleading logarithms whose effect reduces to an overall factor 1 + 0(cq). 
Thus, there is an uncertainty of order 0(eq) in the normalization of the end-point structure 
which must (and can) be eliminated by further calculations. 

The final predictions for the spectral distributions, ready for comparison with exper- 
imental data, are obtained by superimposing the perturbative corrections with the non- 
perturbative effects I have discussed previously. The primordial shape function is smeared, 
of course, (in a well-defined way): its basic feature - a peak in the end-point domain - stays 
intact. 



12 What is the heavy quark mass? 

So far I have deliberately avoided any explicit definition of the heavy quark mass itiq, as well 
as the related parameter A, although both are the key parameters of the theory, appearing 
in virtually every expression. If you start thinking on this issue you will realize that the 
question is not that simple as it might seem at first sight and as it sometimes presented 
in the literature. The distinction between different definitions of mass is obviously a pre- 
asymptotic effect since the difference is of order Aqcd, a A-qcd/^q effect. Since we are 
aimed, however, at a systematic theory of pre-asymptotic effects ensuring accuracy 1/rng 
and higher we have to address the question of mass from scientific positions |35| . 

Isolated quarks do not exist in nature; therefore, the notion of an on-shell quark mass 
is meaningless. Heavy quarks are always surrounded by a light cloud. How much of it can 
one peel off? 

The quantity most commonly used to parametrize tuq is the so-called pole mass fn^ le . 
In perturbation theory the quark Green function is well defined in each (finite) order; the 
pole mass is, by definition, the position of the pole of the quark propagator to a given 
(finite) order in a s . Thus, what people actually do is equivalent to separating out the 
non-perturbative component of the cloud from its perturbative (Coulomb) part. The per- 
turbative component is supposed to be included in itiq while the non-perturbative one to 
A. 

In purely perturbative calculations such a definition is quite acceptable, as well as many 
others. Moreover, it is even more convenient than others since the gauge invariance of rn P Q le 
is explicit. This feature made it very popular since the mid-seventies (see e.g. Ref. 36). 
Important results of the perturbative calculations, such as the total semileptonic widths, 
are routinely expressed in terms of the pole mass. 

What is acceptable in PQCD becomes absolutely unacceptable in the theory of non- 
perturbative effects. As a matter of fact, it is impossible to define fn P Q le to the accuracy 
better than AqcdI^q [35]. The discrimination "perturbative versus non-perturbative" 



contradicts the spirit and the letter of the Wilsonian operator product expansion and can 
not be carried out in a systematic way. The right thing to do is to separate out the low- 
frequency component of the cloud {uj < fi) from the high-frequency one (u> > fi), and include 
the high-frequency component in the definition of the quark mass. Then the low- frequency 
component will reside in A. In other words, any scientific approach must deal with mg(/i), 
the running mass introduced at a distance [i below the would-be pole. As usual, the choice 
of the parameter \i is dictated by two opposite requirements: on one hand, to come as close 
to the soft non-perturbative domain as possible we want /j, to be as small as possible; on 
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the other hand, to develop a consistent theory we must ensure that a s {fj)/iT <C 1. The fact 
that both requirements can be met simultaneously is a highly non-trivial feature of QCD, 
established so far only at an empirical level. Practically \i = several units xAqcd- 

The fact that it is impossible to define m P q le to the accuracy better than Aqcd/ttiq is 
seen in many different ways. One of the simplest arguments comes from the consideration 
of the divergence of perturbation theory in high orders due to the infrared renormalons 
[37-39]. Of course, nobody claims that the presence of the infrared renormalons is the 
only obstruction to introducing m P Q le ', the only uncontrollable non-perturbative piece. But 

their occurrence is quite sufficient to kill any consistent calculation dealing with m P Q le and 
aimed at accuracy AqcdI^Q or better. The renormalon contribution gives an idea of 
the theoretical uncertainty which can not be eliminated, in principle, without introducing 
m Q (n). 

Let us examine the perturbative series for the heavy quark mass. The renormalon 
divergence originates from the "bubble" insertions into, say, the one loop graph (Fig. 11). 
A standard expression for the one-loop renormalization can be written as 

87ra s f d 3 k 1 



. , 87ra s f ark 1 



where (i is the normalization point of the running mass. Summing up all the bubbles we 
effectively substitute the bare coupling constant by the running one in the integrand, 

pole i v 8tt f d 3 k a s (P) 
m Q -mQ(n) = — r . (38) 

v 3 J\k\<v k 2 

Now, the running gauge coupling is given by 

a s (P) = ^ _ = a s (S) t ( * " 

l-(&a a (M 2 )/47r)M M 7fc2) ^ ' tA ^ W 

where b is the first coefficient in the Gell-Mann-Low function. Substituting this expansion 
in eq. (38) and doing the k integration we arrive at a factorially divergent series 



r--»«(f>=^£»>(^) . (3 9 ) 



This series is not Borel-summable; truncating it at an optimal value of n, n$ ~ 2-7r/(&a s (//)), 
we find that the intrinsic uncertainty in m P Q le — raq(fi) is 

A ( m Q le ~ m Q (n)) ~ ^^exp{-— ^-y} ~ ^QCD- (40) 

(For a related discussion of the infrared renormalons in connection with HQET, with similar 
conclusions, see Ref. 40). 

Eq. (38) counts the energy of the would-be Coulomb tail around a static color charge. 
The tail actually does not exist in the non-abelian theory; this is the reason explaining the 
intrinsic uncertainty (40). 

The fact that there is no scientific way to define m P Q le to better accuracy does not mean, 
of course, that one can not construct a systematic expansion in I/thq. This expansion, 
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however, must be based on HQET/OPE and explicitly include the normalization point ji. 
Then the expansion parameter is actually mq(ii). 

Although raq(/j) is not directly measurable it is related, nevertheless, to observable 
quantities. Let us elucidate the last point in more detail. 

Assume that the structure functions in the semileptonic decay B —* X c lu are measured 
separately as functions of qo and |g|; we are specifically interested in w\. (I will not bother 
you with precise definitions, w\ is just one of five possible structure functions, see e.g. Ref. 

Assume furthermore that the measurement can be made at a fixed and small value 
of |g|, i.e. <C Mp. Then the plot of w\ will schematically look as depicted on Fig. 
12. Qualitatively it is similar to the plot of Fig. 8 where E is now replaced by c/o> a long 
radiative tail to the left of the shoulder was omitted on Fig. 8. This tail has nothing to do 
with the primordial motion, it is due to the hard gluon emission. The height of elastic delta 
function is 1, in certain units; a few conspicuous excitations in the shoulder have heights 
0{v 2 ) where v = \q\/Mf) while the height of the radiative tail is further suppressed by 

a S (qOmax - <7o)/tT- 



Using our approach we can formally show [41| that the integral 



j (qomax ~ qo)dq wi = A (41) 



L D 



modulo corrections of higher order in v and in Aqcd- The upper limit of integration is 
Qomax, the lower limit is formally zero. 

In these terms the standard definition of a "universal constant" A, and the correspond- 
ing pole mass, mf^ le = Mb — A, might be formulated as follows: (i) take the radiative 
perturbative tail to the left of the shoulder and extrapolate it all the way to the point 
qo = qomax ; (h) subtract the result from the measured structure function wi ; (iii) integrate 
the difference over dqo with the weight function (qo ma x — Q'o)/2vr. The elastic peak drops 
out and the remaining integral is equal to (v 2 /2)A. 

When I rephrase the standard program that way its absurdity, from purely theoretical 
point of view, becomes evident. How can one define what the perturbative tail is at (^omai — 
qo) ~ Aqcd even in the first order in a s , to say nothing about high orders? Meanwhile, this 
is exactly what is implied in the pole mass: one includes in the pole mass, by definition, all 
perturbative corrections coming from all virtual momenta, up to the would-be mass shell 
(i.e. excludes them from A). 

What is the correct procedure? We must introduce a normalization point \x in such 
a way that all frequencies smaller than \i can be considered as "soft" or inherent to the 
bound state wave function; at the same time a s (/i) is still a well-defined notion and the 
perturbative expansion in a s (n)/ir makes sense. We then draw a line at go = qomax — ^ (Fig. 
12). The integral (41) taken in the limits from qomax — ^ to qomax represents (w 2 /2)A(/x) 
modulo corrections of higher order in v and in Aqcd- The running mass is then defined 
as m&(/i) = Mb — A(/x). Practically, if the radiative tail in w\ is numerically suppressed 
compared to the height of the shoulder, the \i dependence of A may be rather weak. 



13 The onset of duality. Divergence of OPE 

In this part of my talk I will leave particular applications and address a general issue, the 
foundation of the theoretical approach developed. The theory of preasymptotic effects in 
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weak decays of heavy flavors is based on the Wilsonian operator product expansion. More 
exactly, we exploit a generalization of OPE relevant to the Minkowski kinematics. To justify 
the procedure we keep in mind a possibility of analytic continuation to the euclidean domain 
in some momenta. In the semileptonic decays we can analytically continue in the lepton pair 
momentum and consider the amplitudes of interest off all cuts, as it is done, for instance, 
in Ref. 10. In the general case it is always possible to introduce an auxiliary (complex) 
momentum flowing in (or out) the vertex generating the weak decay. The operator product 
expansion is unambiguously defined in the euclidean domain since here the separation of the 
large and short distance contributions is well-defined. Predictions in the Minkowski domain 
are then obtained through dispersion relations. The simplest prediction of this type is the 
calculation of the total e + e~ annihilation cross section at high energies. 

The coefficient functions in the original (euclidean) version of OPE are determined 
by the short distance dynamics. From the remark above it must be clear that the very 
same statement is valid for the physical observables in the Minkowski region provided they 
are fully inclusive and the predictions are understood in an integral sense, smeared over 
some energy range. A point-by-point prediction becomes possible only at sufficiently high 
energies, when the observables under consideration are already smooth and additional "by 
hand" smearing is not needed. The theory of the preasymptotic effects I have discussed 
above assumes that we are already inside this duality domain. The question of where the 
boundary of the duality domain lies in each particular process is one of the major practical 
questions in the non-perturbative QCD. 

To make the point more transparent let us turn to the e + e~ annihilation cross section 
where the kinematical conditions are much simpler than in the heavy flavor decays. If q is the 
total momentum of the pair, q 2 > corresponds to the physical (Minkowski) region while 
negative values of q 2 (i.e. positive Q 2 = —q 2 ) is the euclidean domain. The operator product 
expansion in the euclidean domain is the expansion of the current correlation function in 
the powers of a s (Q) and Aq CD /Q 2 modulo logarithms of the same parameter. As well- 
known J37| , |39| the perturbative series is factorially divergent and non-Borel-summable. 
The standard perturbation theory has to be amended. 

It was amended |42] in the spirit of OPE - by means of introduction of the vacuum 
condensates, 

{vac\O n \vac) ~ \i n exp 



, , logs, (42) 

where d n stands for the dimension of the operator O n . This rather straightforward extension 
cures the divergence of the perturbative series; the exponential (in l/a s ) terms appear 
explicitly. They are obviously non-analytic in a s . These terms result in a complicated 
analytic structure in the complex a s plane, with an infinite set of cuts [37|. 

The expansion for the current correlation function at Q 2 > is translated in the pre- 
diction for the total cross section, o~(s), where s = q 2 . In the limit s — ► oo this prediction 
is simple: an appropriately normalized ratio R = Const (sex) is equal to unity (I assume for 
simplicity that the only quarks are the massless u and d). Of course, we are interested in 
deviations from this unity. 

Every given term in the euclidean expansion results in a correction in R{s) of a generic 
form 

~]a n (A 2 QCD /s) d ^ 2 x logs. (43) 



(Pure powers of 1/Q 2 yield delta function of s and its derivatives invisible at s > Aq CD .) 
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We have calculated a variety of such terms in the case of the heavy flavor decays where the 
role of s is played by nig. It is obvious that at s > A-qcd t° anv finite order in eq. (43) 
the prediction for R(s) is a smooth function of s. 

The question is where the structures in R(s) come from. The fact that the structures 
are there is clearly seen, say, in the multicolor QCD f43]| , with N c — > do. In this limit R(s) 
actually represents a comb of infinitely narrow peaks: the lowest one is due to the p meson 
and the rest are due to its radial excitations. Even if the number of colors is three the first 
excitation shows up in R(s) as a noticeable structure although numerically mf,i/A-QCD ^ 1- 

The answer to the question above is in the behavior of the power expansion (43) in high 
orders [Q] (or, better to say, in the behavior of the analogous euclidean expansion). The 
power series per se turns out to be asymptotic, it has factorially divergent coefficients. The 
reflection of this asymptotic nature of the expansion is the occurrence of the exponential 
terms of the type 

eM-CQ 2 /A 2 QCD ) 

in the euclidean domain. Their analytic continuation to the Minkowski region provides a 
structure in R(s). They also set the rate of approach to the duality domain and determine 
the accuracy of the duality relations. 

The factorial divergence of the power series may be called the divergence of the second 
generation. Accumulation of the 't Hooft singularities (43) produces a more sophisticated 
singularity Q 



exp 



(44) 



Unlike the 't Hooft singularities which, in principle, can be revealed by analyzing multi-loop 
Feynman graphs, the second-generation singularity (44) obviously remains undetectable 
diagrammatically (in the classical understanding of diagrammar). 

The assertions above can be demonstrated in a few different ways. The simplest proof 
I know is as follows. Let us consider heavy-to-light quark currents 

J s = Qq, Jp = 075<? ; 

Jl = \{J S + Jp)= Q^(l+75)9, J2 = \(ji-J P )=q\{l + ls)Q. (45) 

Later on we will assume that uiq —* 00 and m q — > 0. Note that J2 is not hermitean 
conjugate to J\\ rather J2 is a chiral partner to j\. Therefore, the correlation function 

n = ij e fc dx(wac|r{Ja(x), J 2 (0)}\vac) (46) 

vanishes in perturbation theory if m q is put to zero. Thus, there is no uninteresting back- 
ground, a feature which should be welcome, of course. 

Formally we can write the correlation function (46) as a trace, 

n(f)-m {^+ ^ 2 _ m r; (i/2M; (V + k? -rl + m)aa \ (47> 

where V is the momentum operator and the trace operation is defined with respect to this 
operator. Simplifying a little bit I will say that II is the product of two quark Green's 
functions in the background gluon field; the averaging over this field is implied (Fig. 13a). 
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The (1 + 75) projector in the currents J\ 2 ensures the "softness" of V: only very small 
eigenvalues of V 2 + (i/2)aG contribute since otherwise the explicit factor m q in eq. (47) 
annihilates the result in the limit m q — ► (compare with Ref. [pp. 

Let us focus on the soft gluon fields discarding the effect of the hard gluon exchanges 
which is well- understood. First of all, we should choose our reference point, the external 
momentum k, in a most advantageous way. If mg — > 00 it is clear that the optimal reference 
point lies "slightly" below threshold, 

h = rriQ - e, k = 0, 

where e scales with Aqcd, not rriQ. 

Taking now the limit toq — ► 00, assuming the "softness" of V (i.e. V ~ Aqcd] this 
assumption is clearly not valid for graphs with the hard gluons which are disregarded 
anyway) and neglecting all terms suppressed by powers of 1/toq we find that 

The momentum operator V is a perfect analog of the operator ir from the previous sections. 
As a matter of fact, to be fully consistent, I should have replaced it by 7r; and I will from 
now on. If e 3> Aqcd one can expand eq. (48) in 1/e. To write the expansion in a compact 
form it is convenient to introduce the vacuum expectation value of qq where q is the massless 
quark field. Then, according to 



(vac\qq\vac) = lim (-i)Tr < — ^ k — 9 ,. lrk , _ > . (49) 

x 11 ' m 9 ^o v ; [tt 2 - m\ + {i/2)<jG J v ; 

Let us agree that in this section the angle brackets will denote the vacuum averaging. Then 



m = h 



(QQ) + -3(0*09) + ^(Q^oQ) 



+ ... 



(50) 



The term with 75 in eq. (48) drops out since the vacuum state in QCD is P-even. What is 
more important - a key point in our analysis - is the disappearance of the odd powers of 
ttq. The vacuum expectation values with the odd powers of ttq all vanish due to the Lorentz 
invariance. Another way to see the absence of the odd powers of ttq is to return to eq. 
(48) and to observe that it must stay intact under the change of the sign of the momentum 
operator, an obvious property of the trace operation. Certainly, one could have written the 
operator product expansion (50) directly by cutting the light-quark line on Fig. 13a (see 
Fig. 136). 

A question that immediately comes to one's mind is how eq. (50) can possibly be correct. 
Indeed, let us assume that the expansion (50) is convergent. Then it defines a function which 
is odd under the reflection e — ► — e. By no means one can allow the correlation function 
11(e) be odd! 

At positive e we are below the cut, and 11(e) is analytic in e. At negative e, however, we 
sit right on the cut generated by the intermediate physical states produced by the currents 
Ji 5 2, and the correlation function 11(e) develops an imaginary part, a discontinuity across 
the cut. Qualitatively we have a pretty good idea of how ImLT looks like. I will dwell on 
this point later and now just state that on physical grounds the behavior of n(e) at positive 
and negative values of e is absolutely different. 
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What is the solution of this apparent paradox? 

Looking at the definition of the function LT and the currents (45) we observe that LT is 
actually a difference between the scalar and pseudoscalar channels, 



noc (J S ,4) 



(j P ,4 



i r c 
n ( e ) = " / 

7T J Ei 



(51) 



Therefore, at first sight one might try to say that the spectral densities in these channels 
are identical and cancel each other in the difference so that n(e) has no imaginary part 
whatsoever. The expansion (50) itself tells us, however, that this conjecture is wrong. A 
brief reflection leads one to an oscillating imaginary part schematically depicted on Fig. 14. 
The oscillations are needed in order to kill all even terms in 1/e in the expansion of the 
dispersion representation, 

lmU(E) 

IE e + E 

The spectral density on Fig. 14 refers to the real QCD, with three colors (I remind that 
it is the difference between the spectral densities in the scalar and pseudoscalar channels 
and, therefore, it need not be positive). In the multicolor QCD, with N c — > oo, we would 
have, instead, two combs of infinitely narrow peaks sitting, back-to back, on top of each 
other. Let us defer for a while the discussion of how the spectral density evolves from 
N c = 3 to N c = oo. 

Returning to the expansion (50), the only logical possibility is to say that the coefficients 
of the 1/e expansion are actually factorially divergent in high orders, so that the expan- 
sion is asymptotic. Then a function defined by this expansion (plus the correct analytical 
properties, see eq. (51)) can well be essentially different at positive and negative values of 
e, although superficially it produces an impression of an odd function of e. (Examples of 
such functions are very well known in mathematics, and I will present one shortly.) In other 
words, the matrix elements (q(7ro) 2n q) at large n behave as 



-i) 



l (qq}(^QCD) n C 2n (2n)\, 



(52) 



where C is a numerical constant. The value of this constant is crucial since it determines 
the onset of the duality regime. 

The appearance of the factor (— l) n can be explained as follows. The vacuum expectation 
value (52) written as an euclidean integral has the form 



(q{vo? n q) = - lim (-l)"Tr 

m„— >0 



2d 



m„ 



it 2 + m 2 q + (i/2)cxG 



(53) 



where the conventions concerning the euclidean 7 and a matrices are borrowed from [fh|. 
The subscript E marks the euclidean quantities. The expression in the square brackets is 
well-defined and is positive-definite. 

Upon reflection it is not difficult to understand that the factorial growth of the conden- 
sates (52) with n is quite natural. Combining eqs. (50) and (52) on one hand with the 
dispersion relation (51), on the other, we conclude that the absolute value of Imll at large 
E is exponential in E, 



\ImU\ ~ exp(-E/(CA QC D)) = exp 



_]_j2Kb-ya s (E)) 



C 



(54) 



25 



The truncated operator product expansion (50) would predict Imll = at large E. Thus, 
the exponential in eq. (54) gives an estimate of the accuracy of the duality relations. The 
number C does indeed play a key role since the duality sets in at E S> CAqcd- 

Speaking in terms of the analytical structure in the a s plane we deal here with a so- 
phisticated singularity of the type exponential of the exponential of l/a s \ 

I have no time to go into further details. Frankly speaking, quantitatively, not so many 
details are worked out. Let me make only two remarks. 

If the number of colors increases from 3 to infinity the spectral density Imll at large E 
experiences a dramatic evolution. A smooth function depicted on Fig. 14 is converted into 
two combs of delta functions - scalar resonances have positive residues while pseudoscalars 
appear in Imll with negative residues; all residues are 0(1). (Certainly, averaging over 
many resonances produces Imll ~ 0.) This means that the terms in the operator product 
expansion (50) which are subleading in 1 JN C must be more singular in n, i.e. their divergence 
in high orders of OPE must be stronger. This aspect is not new, though; the very same 



situation takes place in the ordinary perturbation theory [47|. 

To wind up the topic let me give an example of a function whose properties are very 
similar to those we found for the correlation function (46). If ip(z) is the logarithmic 
derivative of the T function a schematic model for II is as follows: 

K-model = ^{e + ~) ~ He - -) - e _ . (55) 



From the well-known properties of the tp functions 4£] we infer that H mo dei is a sum of 

(56) 



simple poles. All poles are situated at real negative values of e, viz. 

13 5 7 



4' 4' 4' 4' 

and all residues are ±1. The signs alternate: the first pole has negative residue, the second 
positive and so on. In other words, the analytic properties of H m odei are perfect - analytic 
everywhere except the points (56) corresponding to infinitely narrow resonances. 

On the other hand, at large positive e our model correlation function is expandable [] in 
1/e and, moreover, the expansion looks like an odd function, 

1 °° CL 

n=0 ^ ' 

where a n are calculable coefficients. At large n they are related to the Bernoulli numbers 
B n . As well-known, the latter grow as B n ~ (2n)! (see |4^], page 23), and the coefficients 
a n grow in the same manner. Thus we see that H mo dei is a perfect example, it is exactly 
what we want to get in QCD dynamically. 

14 Conclusions 

Kolya Uraltsev recently reminded me about a remark with which Bjorken concluded his 



talk at Les Rencontre de la Valle d'Aoste in 1990 |29j. Bj said then "... within a year or two 



it is quite possible that the language used in describing heavy quark decay phenomenology 

As a matter of fact, the asymptotic expansion (57) is valid everywhere in the complex plane except a 
small angle near the negative real semiaxis, see sects. 3.12 and 3.13 in |48|. 
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will shift away from comparison of data with Model A or Model T, and instead be phrased 
in a language which deals with the importance of a correction of Type X or Type Y...". He 
was slightly inaccurate in two respects: elements of the QCD-based theory of preasymptotic 
effects had existed before 1990, and it took about four years - not one or two - to complete 
the first stage of the theory. In all the rest his anticipation came true - today any serious 
work related to the heavy flavor decays has to start from the 1 / tuq expansion and analysis 
of different corrections. 

There is still a lot of work to be done. On the practical side, the 1/mg expansion has 
to be built at least up to the level of 1/mg in all problems of interest. We must create a 
catalog of the Hq expectation values of relevant operators, much in the same way as the 
vacuum condensates were classified and estimated within the QCD sum rule method. The 
issue of the gluon radiative corrections has to be exhaustively worked out. 

On the theoretical side, the most important question, the heart of the OPE/HQET 
theory, is understanding the quark-hadron duality: when it sets in, what deviations one 
may expect in any given process, etc. These questions are intimately related to the general 
structure of QCD at large distances and the color confinement mechanism. I do not ex- 
clude that attempts to answer these questions emerged from the needs of the applied QCD 
will help to QCD fundamentalists - those who still hope that an analytic solution of the 
(continuos) four-dimensional quantum chromo dynamics is possible. 
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Figure Captions 



Fig. 1. 

Different inclusive transitions of the heavy quark in the quark language. The closed circle 
denotes an effective electroweak vertex. Thick solid line denotes the decaying heavy quark. 

Fig. 2. 
Nuclear (5 decay. 

Fig. 3. 

Graphic representation of the transition operator in the inclusive semileptonic decay. The 
region inside the dashed line is governed by short-distance physics. 

Fig. 4 

The diagrams determining the coefficients of the operators QQ and Oq in the transition 
operator of Fig. 3. 

Fig. 5. 

The photon energy spectrum in the Q — > 57 transition (the parton-model approximation). 
Due to the two-body kinematics the photon line is monochromatic. The final quark is 
assumed to be massless. 

Fig. 6. 

A realistic spectrum in the inclusive hadronic decay Hq — ► X q j. The kinematical boundary 
is shifted to the right of the parton-model line by A/2. The final quark q is assumed to be 
massless. 

Fig. 7. 

Graphic representation of the transition operator in the problem of the end-point spectrum 
(the Born approximation). All quark lines are in the background (soft) gluon field. 

Fig. 8. 

Evolution of the spectrum of Fig. 6 as the mass of the final quark increases (the schematic 
plot refers to m q ~ mg/2). The effects of the hard gluon bremsstrahlung are not included. 

Fig. 9. 

The photon spectrum in the SV limit ( mg — m q <C iTiQ,q)- The effects of the hard gluon 
bremsstrahlung are not included. 

Fig. 10. 

The Roman model for the distribution function relevant to b — > u. The Roman ansatz has 
two fit parameters, the Fermi momentum pp and the spectator mass m sp . The parameter 
p is defined as p = (m sp /pF) 2 . 

Fig. 11. 

Infrared renormalon in the pole mass. 
Fig. 12. 

A schematic plot of the structure function w\ in the inclusive semileptonic decay B — ► X c lv 
in the limit \q\/Mu <C 1. A long perturbative tail due to the hard gluon shake-off is depicted 
to the left of the primordial domain. The dashed line shows an extrapolation of this tail in 



31 



the region where q 0ma x - <7o ~ A. 
Fig. 13. 

Diagrammatic representation for the correlation function (46). (a) The product of the quark 
Green functions in the soft gluon field; (b) Building the operator product expansion: the 
light quark line is soft in the limit m q — ► and can be cut. 

Fig. 14. 

A schematic plot of the spectral density for the correlation function (46). 
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